Optimal estimation of quantum dynamics 
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We construct the optimal strategy for the estimation of an 
unknown unitary transformation U G SU(d). This includes, 
in addition to a convenient measurement on a probe system, 
finding which is the best initial state on which U is to act. 
When U G 517(2), such an optimal strategy can be applied to 
simultaneously estimate both the direction and the strength 
of a magnetic field, and shows how to use a spin 1 /2 particle to 
transmit information about a whole coordinate system instead 
of only a direction in space. 
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Consider an experimental device T> that implements an 
unknown unitary operation U € SU(d). A probe subsys- 
tem A, which can be entangled with a second subsystem 
B, is introduced in V and analyzed at its releasing. Sup- 
pose that arbitrary manipulation is allowed on the global 
composite system both at the preparation and analysis 
stages, while V is regarded as a black box. This paper 
addresses the question: "Which is the best way of esti- 
mating the operation UT' 

The optimal estimation of the state of a quantum sys- 
tem has received a lot of attention in recent years ■ 
A situation repeatedly considered in the literature is that 
of a spin 1/2 system prepared in an unknown pure state 
€ C 2 . By means of an optimal measurement on the 
system, the maximal amount of information about \ip) is 
retrieved. Here we focus, instead, on the estimation of 
the dynamics of a quantum system (see also Q ) . This is 
done by analyzing, again through an adequate measure- 
ment, the changes that the initial state \ipo) S C d <g> C d 
of the system undergoes under the unknown evolution, 
U G SU(d). But contrary to what happens in state es- 
timation, where only optimal measurements need to be 
constructed, the optimal estimation of transformations 
requires a double maximization: first, and most novel, we 
need to find the state \ipo) of the composite system that 
best captures the information of the transformation (uni- 
tary evolution U); and second, a measuring strategy that 
optimally retrieves such information from U <g> Ib\i/}q), 
where / stands for the identity operator. 

Not surprisingly, the optimal estimation of quantum 
transformations — necessarily based on the possibility of 
encoding them on, and analyzing them from, a quantum 
system — is closely related to the capacity of quantum sys- 
tems to carry information. Our results also give insight 



into the role entanglement plays at enhancing the capa- 
bilities of a quantum channel: it turns out that unitary 
transformations are optimally encoded in the quantum 
correlations between the two subsystems, A and B, and 
that, for instance, information about a whole coordinate 
system {e x , e y ,e z } can be transmited by sending only one 
spin 1/2 system, provided that an ebit of entanglement 
between the sender and the receiver is also available. The 
simultaneous determination of both the direction and the 
strength of a magnetic field, the tuning of a quantum 
channel and the limits to espionage in a two-party pro- 
tocol are other issues that can be addressed with the 
optimal scheme for the estimation of unitary operations, 
as we shall discuss. 

It is easy to come up with strategies that determine 
U with an arbitrary accuracy provided that the black- 
box device V can be used without restrictions. Here we 
are interested in the opposite situation, namely when D 
is used to perform the transformation U only a reduced 
number of times N. We will first present an exhaustive 
analysis, comprising the optimal initial state |^o) and the 
optimal measurement, for the case when T> can only be 
used once, N = 1. For the general N case, and assuming 
that T> performs the transformations in the form U® N , 
we will derive the optimal initial state of the system, and 
report the optimal POVM for N = 2, U € SU{2). 

We start by shortly reviewing some of the elements in- 
volved in quantum estimation strategies. First, a prior 
probability distribution f(U) uniform with respect to 
the Haar measure || expresses the fact that nothing is 
known about U before resorting to T>, except that it cor- 
responds to a unitary evolution. Second, once the device 
T> has performed U on the probe A, a positive operator- 
valued measurement (POVM) on A and the (possibly) 
entangled system B will extract the information about 
U. Such POVM is a set {G r } of positive operators satis- 
fying J2 r G r = Iab- And third, we need a notion of how 
efficient a particular strategy — that is, an initial probe 
state |V>o) and a POVM {G r } — is, so that we can search 
for the best one. There are several ways of evaluating the 
strategies, and the optimal solution may depend on the 
particular election we make. One of the main results of 
this paper is to present the optimal probe state \ipo) and 
to show that it is the same for a large class of figures of 
merits. Nevertheless, in order to optimize the POVM we 
will consider a specific, fidelity-guided figure of merits, in 
which the outcome r of the POVM, corresponding to the 
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operator G r , is followed by a guess U r for the unknown 
U. We have chosen the function 



F(U, U r ) 
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tr(C/C/t)| 2 (1) 



to evaluate the guess U r . It quantifies, on average over 
all states \4>), how well U r compares to U when trans- 
forming \(f>), since it averages the overlap between U r \(j)} 
and U\4>) . Bellow we will give another interpretation to 
this fidelity, whose average over outcomes and unknown 
operations reads 

F = V/ f(U)dUP r (U)F{U,U r ), (2) 

r JSU(d) 

where P r {U) is the probability that the POVM produces 
the outcome r when the device T> has implemented the 
operation U . 

Let us suppose, then, that V is to be used only once. 
Lemma 1 presents the optimal initial state of the probe 
for this case. It only assumes a covariantly averaged fig- 
ure of merits as in (Q), but where F(U, U r ) is any function 
h(UU}) depending on U and U r through UUj. Notice 
that only pure states need to be considered for the probe 
system, due to the linearity of P r {U) in the initial state 
(see eq. (Q)). Therefore we take, without loss of gen- 
erality, a composite probe AB, where A is the d— level 
system on which U will be performed and B is a second 
d— level system, possibly entangled with A. 

Lemma 1: The optimal initial state for estimating U 
after a single performance can be chosen to be a maxi- 
mally entangled state, such as 
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1>A1>B> 
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The reason is that, as we next show, the state U (8 Is \&) 
can be transformed, independently of U, into any other 
state U (8 Ib\^o) — actually, to an equally efficient state, 
see below — by just manipulating system B. 

Proof: Let us consider the Schmidt decomposition of 
the most general initial state |V>o) = Ylj—i Aj > 

Ai+i > 0, J2i^H = 1- We first show that the Schmidt 
basis {\(J-iVi)} is irrelevant as far as the average fidelity 

h = ^tx(G r J dUU®I B \ip ){i> \tf®I B h(UUl) 

(4) 

is concerned (here ti (G r U <g> I B \'4'o){'4'o\U' < <8 Ib) is the 
probability P r (U)). This is so because for any X and 
Y € SU(d), the state Xa <8 Yb\iPq) leads to the same 
maximal h, as can be seen by noting that: (i) any uni- 
tary transformation Y in the local basis of B can be re- 
absorbed in the POVM elements G r , whereas (ii) if we 



prepare AB in state X ®I B \ip ) instead of |-0o) •> t nen the 
shift U —> UX in the integration variables U of eq. (Q), 
simultaneous to a shift U r — > U r X for the guesses leads 
again to the same /i, as a consequence of the isotropy of 
f(U). Therefore we can take 
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I A ®M\$), 



(5) 



where M is a diagonal operator with entries Ma = 
Aj/Ai < 1. Suppose now that the initial state is |$). 
Then T> transforms it into U ® Ib\$)- Let us con- 
sider a covariant POVM Q on B given by operators 
{M Y = (d/trM 2 )2MF}, where Y runs isotropicaUy over 
SU(d) and / dYM Y M Y = I B - It transforms U ® I B \&) 
into (d/trM 2 )iu®MY\<P) = (d/trM 2 )^UY T ® M|$) = 
L/V t C>5/b|V'o) for some known V [here we have used that, 
VV G SU(d), I A ® V|$) = V T ® 7s|$)]. But this is as 
if we would have started with state Y T ® /s|Vo), which 
leads to the same average fidelity as |Vo)- 1=1 

Let us now notice that our particular choice of fi- 
delity, eq. (0), corresponds precisely to the probability 
\(^\0U(g>I B \^}\ 2 that the state U r ®I B \<&) behaves as if 
it were U ® I B |$) . Therefore i r (?7, C/ r ) measures how sim- 
ilar the two operations U and U r are by comparing two 
related states: those that best capture the information 
of both transformations after a single run of D. 

Suppose finally that system A, in the entangled state 
|$) with system B, has already been introduced in the 
black box T), which produces the state U (8 I B \$) — 
denoted by J7|$) from now on — . Which is the best 
POVM that can be performed in order to learn about 
U from this state? We can rewrite the average fidelity of 
Eq. (|) as 



i = ^E tr G r //(f/)^l < f)<*|C/ t |tr(C/C/ r t )| 5 



(6) 



By means of a shift U — > V — U},U in the integration 
variables, each of the integrals inside the trace has the 
form U r fiUJ, where 



d 2 ($\ / /(y)dT/F 02 (|$)($|) 82 y tcs2 |$). (7) 



fi 



Schur's lemma Q states that this last integral is propor- 
tional to the identity in each of the two corresponding 
irreducible representations of SU(d), namely the sym- 
metric and the antisymmetric ones. A careful analysis 
[recalling that each V is acting only on the first half of 
the corresponding |<f>)] and patient simple algebra leads 
to 
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h = d^T \~l^ lA ® Ib + ■ (8) 

Thus |$) is the eigenvector of fi with greatest eigenvalue, 
X m = 2/d 2 , and tr( U^G r U r f! ) < A m trG r in eq. (§). 
Since trG r = cP, the maximal fidelity can be 2/d 2 at 
most. A covariant POVM || with operators and guesses 
given by {W|$X*|W+, W}wesa(<J) reaches F t = 2/d 2 , 
which is consequently the optimal one. 

This result is to be compared with the optimal fidelity 
Fa = 1/d 2 made by blindly proposing a unitary transfor- 
mation, say I (or any other): 

J f(U)dU 1 -^ = (*| J dUf(U)U\m\U^m (9) 

(the last integral is simply 1/d 2 because of the Schur's 
lemma) and also with the separable fidelity = 
(d + 2)/[(d + l)d 2 ], which is the best fidelity that can 
be achieved without entangling A and B, and can be 
computed using eq. (Q) and the fact that a pure state 
of A, say |0), is \/d(0s|<f>). Finally, we note that a fi- 
nite (and thus physical) optimal measurement, actually 
one with the minimal number of outcomes, consists in a 
von Neumann measurement on a basis of d 2 maximally 
entangled states. For instance, on the Bell basis, with 
guesses I,ia x ,ia y and ia z , for the SU(2) case Q|. This 
completes the analysis of N = 1 (|] . 

Let us discuss some applications of the previous results. 
Consider first the group SU(2). Our optimal strategy 
can be readily applied to determine a constant magnetic 
field B = Brh by using the magnetic moment of a spin 
1/2 particle, say an electron. Let Hi nt = j2 ■ B be the 
interaction Hamiltonian, where jl = ii{a x , a y ,a z ) and all 
physical constants have been absorbed in /j,. Then after a 
time T the spin has evolved according to exp(—i/j,BTrh ■ 
<?), and therefore we can identify the direction rh of the 
magnetic field and its intensity B (actually /j,BT). Our 
results show how to optimally extract information about 
B by means of an electron if this interacts once with the 
magnetic field (see also 

In the discussion above the information about the mag- 
netic field B is not contained in the state of the spin 
alone, but in the correlations between this spin and a 
second one. Similarly, if two distant parties, Alice and 
Bob, want to use a recently established d-dimensional 
quantum channel, 

d d 

^Cj|u> — >^2ci\i B ), (10) 

i=l i=l 

but Bob does not know the correspondence between 
states — that is, he ignores the states {|«_b)} — , they can 
benefit from a maximally entangled state |$) in order 
to tune the channel. Indeed, by Alice sending her half 
of |$) down the channel, Bob can estimate the whole 



unknown basis or, equivalently, the transforma- 

tion U = SjKs){^A|, with a fidelity 2/d 2 , which is 
2(d + l)/(d + 2) times greater than the fidelity he could 
have obtained also after a single use of the channel if 
no entanglement would have been available. In a sense, 
this is a general manifestation of how entanglement en- 
hances the capacity of a quantum channel, with tradi- 
tional quantum super-dense coding appearing as a 
particular case, namely when the channel is used to trans- 
mit classical information only. 

Let us further see this in the SU(2) case, by assum- 
ing that a spin 1/2 particle is used as a channel. Here 
an ebit of entanglement allows to transmit, by sending 
a single spin 1/2 particle, information about a whole 
transformation U(h,w) S SU(2) or, equivalently, a ro- 
tation R(n,u>) € 5*0(3). In other words, instead of us- 
ing the spin of the particle to try to establish a com- 
mon direction n in space (that of the one-qubit pure 
state \ip)(il>\ = 1/2(1 + n ■ a)), Alice can now send in- 
formation about a whole coordinate system {e x ,e y ,e z } 
to Bob in order to establish a common reference frame. 
This works as follows. The parties share the state 
|$) = (|1a1b) + |2a2b))A/2, where {\i A )} and {\i B }} are 
given with respect to reference frames of Alice and Bob 
respectively. Each party knows his/her own reference 
frame, but ignores the other one. If Alice sends her half 
of 1$) to Bob, then Bob can estimate the rotation R(n, u>) 
(or corresponding unitary U = \1b)(1a \ + I 2b) {2a\) that 
relates the two coordinate frames. 

Another scenario in which these results are relevant is 
that of two parties that are to collaborate in some task 
but do not trust each other. For instance, Bob needs to 
compute on a given input state 1^) a function (unitary 
U) that Alice's computer can perform, but he ignores U. 
Alice is willing to assist Bob by computing U\ip), but 
without letting him find out which transformation U is. 
In this case Alice knows that Bob can estimate U at most 
with a fidelity 2/d 2 . 

So far we have analyzed a single run of the device T>. 
In practice, one would like to determine U with arbitrary 
precision, and this is only possible if T> is used many 
times. Suppose U is performed twice. A most general 
strategy consists on sequentially introducing two probes, 
Ai and A2 , on T>, but allowing for an arbitrary manipu- 
lation of the proofs in between. We do not know how to 
tackle the problem in its full generality We will suppose 
ad hoc that the device T> takes N probes, A\...An, and 
transforms them according to U® N . This could corre- 
spond, in the SU (2) case, to letting the spin of N elec- 
trons interact with the constant magnetic field B during 
some time interval T. 

The first step towards an optimal strategy for estimat- 
ing U is again to find an optimal initial state IV'o') for 

the N d- level systems A = _4_i 4.iv and N auxiliary d- 

level systems B = B\...Bn, that lemma 2 presents. The 
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I/®-" representation of SU(d) contains (several copies of) 
q inequivalent irreducible representations (IRREPs), la- 
beled by a = l,...,q in what follows. For each a there 
are n a equivalent IRREPs, labeled by a(3, (3 = 1, n a , 
each one having dimension d a . The set {| a/3 k) de- 
notes an orthonormal basis for the IRREP a/3, P a $ = 

J2 k \a(3k}(a(3k\ and P a = Ea=i P <±P- The a P and a P' erators {Q Y = E Q Ea a ap Efc \<l>af]k}(a(3k\Y® N }, where 



EFMajTLlptPc,,^!, where EiQ?Q? = Pa 
and the sum (3-\-i is modulus n a , takes with certainty the 
state U® N \& N ) into U® N \&), which is still maximally 
entangled in each IRREPS aj3 but with different weights 
a-ap/a-a in each IRREP, where E/3( a a/3) 2 = a 2 a . And 
second, a covariant POVM in B given by the set of op- 



IRREPs being equivalent, there exists a unitary II' 



such that U® N \af3k) = W p0 ,U® N \af3'k) for any U and k 

§• 

Lemma 2: The optimal initial state for estimating 



\*») EE £ 



(ii) 



where the value of a a > depends on the figure of merits 
under consideration and where 
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Y^\apk) A \apk) 1 



(12) 



0=1 k=l 



is a maximally entangled state between the subspace of A 
that carries the n a IRREPs a(3 (i.e., between the support 
of P a ) and an equivalent subspace of B. For instance, for 
the N = d — 2 case, the optimal initial state is 

= «4 E I**>aI*>b + v/I^k) J4) B , (13) 

v- 3 fe=l 

where \t k ) e {|00>, (|01) + |10))/%/2, |11)} are the triplet 
states and |s) = (|01) — |10))/v2 is the singlet state. 

Proof: Being a generalization of that of lemma 1, here 
we will only sketch the proof. Notice that any state \iPq) 

Ea=l IV'a); where 



of the probes can be writen as |^>q 



\^a) = EE \^k) A \cj> a p k ) B 



(14) 



(3=1 fc=l 



is the projection P a (g) Ib\^) and |0 Q /3fc) are arbi- 
trary states of £?. Since [/^^ does not mix IRREPs, 
we can perform a global unitary transformation Vab 
that commutes with U® N and such that we achieve 
{4> a f3k\4> a 'i3'k') = S^a'Sp^'C^, , that is, the supports of 
Pafi ® ^s|^> on B for different IRREPs a(3 and a'/3' 
are orthogonal. For instance, in the d = N = 2 case, 
where \ipf) = £ fe \t k ) A \<t>k) B and IV'f > = \s)\<t>), we can 
take, without loss of generality, (4>\<fii) = 0. We will now 
show that |$ ) can be transformed into a state as effi- 
cient as \ipo) as far as the fidelity 
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G r / dUU m \i>$){^\tf® N h(UU}.) 



is concerned. This is made in two steps. First, 
the POVM in A defined in each a by {Qf = 



(E 



c a/3 l 



_1 / 2 , will produce, 



the state (UY T )f N \^). 



J dYQ\Q Y = Ib and a a p = 
when applied on Uf N \<f>') 

This state corresponds to starting with Y a ^ n \^Jq), 
which leads to the same h as \tpo) (see lemma 1). □ 

For N = d = 2, and by using the techniques de- 
veloped in this paper, we have found that the optimal 
fidelity is F 2 = (3 + y/E)/8 » 0.6545, which corre- 
sponds to the initial state |$ 2 ) of eq. ( p"3| ) with a 2 = 
(5 + V&)/10 and to a covariant POVM and guesses given 
by {W® 2 \$ 2 a ,)($l,\Wi® 2 ,W}, a' 2 =9/10. 

In conclusion, in this letter we have studied the optimal 
estimation of an unknown unitary operation, U € SU(d), 
when this transformation can be performed a reduced 
number of times, N. For any N the best initial state has 
been essentially found for a large class of figures of merits. 
In the case of the fidelity defined in ([!]), its optimal value 
and the measurement that attains it are given for any 
dimension when JV" = 1, and for d = 2 when N = 2. 
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